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The class of gravito-electric, algebraically general, rotating 'silent' dust space-times is studied. 
The main invariant properties are deduced. The number io of functionally independent zero-order 
Riemann invariants satisfies 1 < io < 2. The solutions may be subdivided accordingly, and special 
attention is given to the subclass to = 1. Whereas there are no A-terms comprised in the class, the 
limit for vanishing vorticity leads to two previously derived irrotational dust families with A > 0, 
and the shear-free limit is the Godel universe. 

PACS numbers: 04.20.-q, 04.20.Jb, 04.40.Nr 



oo 

o 

O 
(N 

1—5 



I 



> 

(N 
(N 

cn 

o 

oo 

o 



X 



1. INTRODUCTION 

Dust space-times in general relativity, characterized as 
perfect fluids with vanishing pressure p = and non- 
vanishing matter density /i 7^ 0, with or without cos- 
mological constant A, have attracted a lot of attention. 
Whereas irrotational dust models form important arenas 
for studying both the late universe ^^ and gravitational 
collapse H , exact rotating dust solutions may serve to de- 
scribe phenomena on a galactic scale. Due to conserva- 
tion of momentum the dust flow lines must be geodesies, 
and the remaining kinematic variables are the expan- 
sion scalar 6 = u°'-a, shear tensor cTab = W(a;&) ~ \Ohab 
(9ab = M(a;fc) being the expansion tensor) and vorticity 
(or rotation) vector w" = \e°'^'^Ub-^c, where u°' is the nor- 
malized dust 4-velocity, hab = Qab + UaUb the projection 
tensor onto the comoving rest space, gab the space-time 
metric and €abc = Vabcdu'^ the spatial projection of the 
space-time permutation tensor rjabcd- 

The Weyl tensor Cabcd, representing the locally free 
gravitational field, is fully determined by its electric part 
Eab and magnetic part Hab w.r.t. u°, defined by 



Eab — Cacbdu'^U , Hab — TT^a 



iC bdu 



(1) 



When Hab = 0, Eab 7^ the dust space-time is called 
gravito-electric, and the Petrov type is necessarily I or D. 
As the gravito-electric tensor is the general relativistic 
generalization of the tidal tensor in Newtonian theory 
(see e.g. 3), such space-times were alternatively termed 
'Newtonian-like' in [5]. 

For Newtonian-like irrotational dust models, the co- 
variant propagation equations for ffab, 0, Eab and fi, 
(equations (|llj|12p and (|14j|15p below, with uj'^ = 0) form 
an autonomous first order system, leading to a set of 
ordinary differential equations when projected onto the 
common eigenframe of (T°b and £'"&. No spatial gradients 
appear, such that each fluid clement evolves as a sepa- 
rate universe, once the constraint equations are satisfled 
by the initial data. Such models were therefore termed 
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'silent' in [6[. The setting looked very appealing towards 
numerical schemes and simulations in astrophysical and 
cosmological context, e.g. for the description of structure 
formation in the universe and the study of the gravita- 
tional instability mechanism in general relativity [7|, |8[ , 
where a clear motivation for taking w" = 0, Hab ~ was 
given in [9]. 

However, in two independent papers [lOJ, lUj, the prop- 
agation of the constraint Hab = along u° was shown to 
give rise to an inflnite chain of integrability conditions 
and corresponding constraints on the above autonomous 
system. These constraints are identically satisfled for 
Petrov type D, but led the authors to conjecture that for 
the algebraically general case only orthogonally spatially 
homogeneous solutions of Bianchi type I (i.e., Saunders' 
cosmological models [l2|) would be allowed. In [ij], how- 
ever, the conjecture was shown to be false for strictly 
positive cosmological constant A, by the explicit con- 
struction of two solution families characterized by the 
presence of a geodesic space-like Weyl principal vector 
fleld. The generalized conjecture that these two families 
exhaust the inhomogeneous Newtonian-like ID models of 
Petrov type / was put forward in [ij| [41[ . 

The idea behind 1+3 covariantly silent models was ex- 
tensively explained and deepened in the introduction of 
[lO|, while in the discussion section of the same paper 
weaker conditions than Hab = tv"" = were indicated for 
establishing the silent property. One of them is to al- 
low for vorticity. This is a natural generalization, since 
silent perfect fluids must have a vanishing spatial gradi- 
ent of pressure [lO| and hence are non-rotating or dust 
(as follows from the Frobenius theorem |16|). However, it 
was questioned at the same time whether this would ap- 
preciably broaden the class of silent solutions, as severe 
restrictions at flrst sight remain. 

On the other hand, important classes of rotating dust 
models have been found by assuming some kind of sym- 
metry, or are algebraically special. Respective examples 
are Winicour's classiflcation 17] of stationary axisym- 
metric models satisfying the circularity condition (see 
e.g. [iMl)) and the general rotating dust solution admit- 
ting time-Uke conformally flat hypersurfaccs with zero 
extrinsic and constant intrinsic curvature as found by 
Stephani [10| and generalized by Barnes for non-zero 



A [2Cf|, which depends on seven free functions of one 
coordinate and which turns out to be gravito-electric 
and of Petrov type D. However, no algebraically gen- 
eral and asymmetric rotating dust solutions have so far 
been found [42] . 

The above serves as a clear motivation to investigate 
the class of Petrov type /, gravito-electric rotating dust 
models, which will be denoted by A, in more detail. In 
this paper the main invariant properties of such models 
are deduced and a natural subdivision based on the num- 
ber ig of functionally independent zero-order Riemann 
invariants is made. The paper is organized as follows. 
Section 2 introduces the 1-1-3 covariant and Weyl princi- 
pal (WP) tetrad settings. The main result is contained 
in section 3, where it is proved that each member of A 
necessarily has A 7^ and possesses a space-like geodesic 
Weyl principal vector field which is parallel to the vor- 
ticity, and where a completed set of algebraic relations 
between the basic scalar invariants (WP tetrad curvature 
components and connection coefficients) is presented. In 
section 4 further invariant properties are deduced and it 
is shown that 1 < to < 2. The subclass corresponding 
to io = 1 is partitioned into two families: the first one 
comprises all non-expanding members of A (for which 
necessarily A < 0), whereas the second one is character- 
ized by equal norms of the vorticity and shear tensors (for 
which necessarily A > 0). Section 5 treats the zero vor- 
ticity, shear-free and Einstein space limit cases. The final 
section summarizes the results in a theorem and briefly 
discusses further features. 



2. MATHEMATICAL SETTING 

We use geometric units SttG — c — 1 and the signa- 
ture (— , -I-, -I-, -I-) for space-time metrics. Round (square) 
brackets denote (anti)symmetrization. By definition, a 
space-time belongs to A if its metric gab is a solution of 
the field equation 

Rab- ^R9ab+J^gab = fJ'UaUt, fl ^ 0, (2) 

and if 



where 



Hab = 0, T2^ - T32 ^0, c^ V 0, 



T2 = &E\E\, T3 = mE\E\E\ 



(3) 



(4) 



are convenient multiples of the quadratic and cubic Weyl 
invariants / and J [18]. 

In the l-t-3 covariant approach [2l|, [24], the tensorial 
quantities uJa, (Jab, 0, Eab and fi play the role of the fun- 
damental dynamical fields, while the parameter A and 
the covariantly constant fields hab and eabc also enter the 
governing equations. A dot denotes covariant derivation 



along M° ('time propagation'). For convenience we de- 
fine, for arbitrary tensors T°'''"cd... and for any natural 
number n, the operator On by 



On [Tf'- • • cd. . . = T'^"- ■ ■ ,d. . . + nOT-'- • • ,rf. . . . (5) 



We further use the streamlined notation of ^23]. The 
spatially projected, symmetric and trace-free part of a 
tensor Sab is denoted by 



1 



S{ab) = ha'^hb'^S(cd) — IzScdh"^ hab — 0. (6) 

The covariant spatial derivative Da, acting on arbitrary 
tensors T°-^-- cd..., and the associated curl and divergence 
(div) operators, acting on one-tensors Va and two-tensors 
Sab, are defined by 

i^eT-^'-crf.. = h%h\ ■ ■ ■ h/h/ ■ ■ ■ hjTP'i-rs.,f, (7) 

div V = DaV", curl Va = eabcD'' V\ (8) 

div Sa = D^'Sab, curl Sab = e.^iaD' Sb)'' ■ (9) 

The Ricci identity for u°- and the second Bianchi iden- 
tity, both incorporating the field equation ^ via the 
substitution 



R. 



'ab = f - + AJ hab + f - - AJ UaUb, 



(10) 



are covariantly split into time propagation and constraint 
equations. For members of A these are: 

• Time propagation equations: 

Ol (9) = -dab'j''' + ^U-^UJa - J + A, (11) 

C| {cr)ab = -CTc{aO-b)^ - '^(a^b) - Eab: (12) 

O^_{L0)a^(TabUj\ (13) 

Oi{E)ab = 3crc(a£^b>'' - ^''(■cdiaEb)'^ - -zdab, (14) 



Oi(/i)-0. 

• Constraint equations 
2 



(15) 



-DaO ~ div (Ta + curlwa. 



(16) 

(17) 
(18) 

(19) 
curl£;„b==0, (20) 

C(l), = -iEabJ" - {a, E\a + [LUJa - 0, (21) 



CUrlo-Qb -I- Di^a'^b) = Hab 

divw = 0, 

div^a - -Dafi = 0, 



where [a,E]a = eabc<^^ dE'^'^ is the one-tensor dual 
to the commutator of a°'b and E°-b- 

Notice that the set of first order equations (|imi5p . aug- 
mented with A = hab = ^abc = 0, indeed forms a 'silent' 
dynamical system, which is due to u°' = cxaXH ab = 0. 



Expressing Hab = and div-ffa = 0, on the other hand, 
respectively yields the two constraints (^0]) and (pij) . 
the repeated time propagation of which gives rise to two 
chains of integrability conditions. In the next section 
attention will be focusscd on the second one, leading to 
a completed set of invariant relations for A under which 
the first chain turns out to be identically satisfied. 

For a generic model contained in A, let B = (9o° = 
u°- , di°- , 82°' , d^"-) denote the WP tetrad, i.e., the essen- 
tially unique orthonormal eigenframe of £'"{, {E12 ~ 
Ei3 = E23 = 0). Below, capital Latin letters A,B,... de- 
note WP tetrad indices and run from to 3, while Greek 
letters a,f3, . . . are {di'^, 82°", 93")-triad indices, run from 
I to 3 and have to be read 'modulo 3' (e.g. aa+ia-i = o'i2 
for a = 3). Einstein's summation convention applies for 
both kinds of indices. In particular we have hap = Saf3, 
with S the Kronecker delta symbol, and by convention 
we take £123 = 1 [231 • The action of a vector field X"' on 
a scalar function / is denoted by Xf. The commutator 
coefRcients 'J'^bc ^^'^ Ricci-rotation coefficients T'^bc of 
B are defined by 



[dB, 8c] = l^BcdA, {dB)^;C = r^BC 



(22) 



As for any rigid frame, the lowered coefficients Tabc = 

gAD^'^BC = ^[AB]C and JABC = 9ADl° BC = lA{BC] 

are biunivocally related by 

lABC = -'2Ta[BC], ^ABC = -1[AB]C + -^IC'AB- (23) 

Within the orthonormal tetrad formalism w.r.t. B, the 
matter density /i and suitable linear combinations of the 
eigenvalues Ea of E"'i, and of the rotation coefficients 
^ABC play the role of invariantly defined basic variables. 
In the present paper we will use 



(24) 
(25) 
(26) 



^a — ro[a-lQ+l]j 
fa+lQ-l = ro(Q-la+l), 

= Pa + iOa+l — r„_iOa-l, (27) 

Xa = i?Q+l— i?Q-l, (28) 

where h\ = — (/i2 -t- /13) and x\ = — (2^2 + 2^3), together 
with [43] 

^a ^ Pq-Iq+io, (29) 

Ua ^ Pa + la-lQ, (30) 

^_a ^ J-aa — la — 1 ^ Ta — 1 a — la \^^) 

fa = Pqq + 1q+1 = 7a + laa+l. (32) 

In terms of the Xa the quadratic and cubic Weyl invari- 
ants read 

T2 = 4:{xl+X2X3+xl), T3 = -4(a::i-a;2)(a;2-a;3)(a;3-a::i) 

(33) 



and the Petrov type / condition T2 — T^ 7^ becomes 
a;iX2a::3 ^ 0. 

The basic equations of the formalism are the commu- 
tator relations, i.e., the first part of (|22p applied to scalar 
functions / (further denoted by cothbc /)j and the pro- 
jected Ricci and second Bianchi equations. Within the 
ON formalism the Ricci equations are further split into: 

1. the first Bianchi equations Raibcd] = 0, equivalent 
to the Jacobi identities 

[diA, [8b, 80]]]"" = 8[ai''bc] + 1^[bc1a]f = 0- (34) 

Here [d[Q, [9^+1, 9„_i]]]° = is the a-component of 
Uni), while [^[1, [92,93]]]° = is (HID; 

2. the a;0-components of (HH) and (fTT]) : 

3. the tetrad components of the field equation ^. 
Here the 00-component is Raychaudhuri's equa- 
tion (|lip and the Oa-component is (fT6|) : the af3- 
components 

8c^ (a/3) — 9(^r a)C + T DC^ (a/3) " T aD^ pC 



= (|+A)5a, 



(35) 



are not covered by the Ricci- identity for w°. 
The formulas 

Va = 8oVa + ea0jn^V\ (36) 

Safl — doSa/S + '^^■ySia^^ Sp) , (37) 

Dcyp= do^Vp-VsT^po., (38) 

DaSfjj — 8aSp^ — 2S'5(^r ^)q, (39) 

and definitions ([Sjn]) relate the covariant differential op- 
erations to directional derivatives. 

3. MAIN RESULT 

Starting from C^^' a ~ one subsequently derives the 
necessary conditions 
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10 



-CTa Eb UJc + -IT-Ea CTb ^c - CTfe E^ LJa^Q 



and 



C(3), = Oz(C«)a 



4 35 

= ■T'^a (Jb'^Ec LOd + -TT^^a E^'^Uc l^d + ^Ea ab'^Uc UJd 

-2at,'E/sd''cOa - 7a/Ed'aa''LUb + 2a/ad'Ej'u:b 

4 5 

--UJ^UJcEa^UJb - i;Ebc^^^''^a - QEaEb'uJc 

\ T7I C TTi b o be I r^ C b 

+ Eb Ec LUa - 6li(Ja (Jb ^c + -^CTb CTc ^a 

5 5 
--eafccw''^'dO-''ew' + -Ej'ebcdto^a'^eto^ 

-t;<TJ'tbcdUJ^E\uj^ +JtbcdE\a^''iJa = 0. 



The validity of these equations has been independently 
checked in an unspecified tetrad approach, hereby using 
the Maple computer algebra package. Note that the ex- 
pansion scalar 9 does not enter the expressions, which is 
a great technical advantage in view of later elimination 
processes. 

In this section we will show that the chain of alge- 
braic integrability conditions generated by further time 
propagation of C'-'^-'a = terminates. We will eventually 
describe its complete solution set both covariantly and in 
terms of WP tetrad invariants. This will be achieved by 
a number of propositions. 

Denote 

Fa = ha+lXa-1 + h^^iXa+l, (40) 

Za = fl{Xa+l - Xa-l) + '^{Xa+l + X^^i). (41) 

Projection of C'-'^'a = w.r.t. the WP tetrad gives 



'a-l-lQ-l 



(42) 



and substituting this in the components of C'-^'' a — 
yields 



u}iX2X3Fi + 2a;2W3a;i-^i = 0, 
UI2X3X1F2 + 2LJ3W1X2Z2 = 0, 
U3X1X2F3 + 2UJ1W2X3Z3 = 0. 



(43) 
(44) 
(45) 



The following lemmas allow to draw quick conclusions 
in later proofs. Lemma [3Tl] especially helps to avoid ex- 
plicit calculations; for its proof we need (fT5|) together 
with the diagonal components of (|14p . namely 



^^ 



Oi{x2) = {x2 + X3 - ^\ h2 + 2:2/13, 
Ci(a;3) = -2:3/12 - [x2 +X3 + -^j hs 



(46) 
(47) 



Lemma 3.1 Suppose 2:2, 2:3 and fj. are constrained by 
two relations 

F{x2,X3, fi) = G{x2,X3,^) = 0, where F and G are 
homogeneous polynomials with integer coefficients and 
without common factors. Then either /i2 = /13 = or 

/i'=T2. 

Proof. As F and G do not have a common factor, 
their resultant w.r.t. 2:3 is non-zero [44] and thus leads to 
at least one irreducible homogeneous polynomial relation 
P{x2,y) = 0. Let n be the total degree of P. As Oi(/i) — 
0, we obtain 0„(P(x2,m)) = ^(2:2, m)Ci(2;2) =0. Now 
■i^ and P are both homogeneous and cannot have a 
common factor, since ^^ has a strictly lower degree than 
P and P is irreducible. Hence ^£-(x2,/i) = P(x2,/i) = 
would lead to (2:2, /i) = (0,0), which is excluded. Thus 
£'1(2:2) = 0, and an analogous reasoning based on the 
resultant of F and G w.r.t. X2 yields 01(2:3) — 0. By 
(|46ll47p this gives a system of two linear and homogeneous 
equations in (/i2:/i3)- Hence either /12 = ^3 = or the 
determinant of the system matrix, computed to be (/^^ — 
r2)/4, vanishes. D 



Lemma 3.2 (a) If two Fp^s vanish at the same time, 
then ft.2 = /13 = 0. 

(b) Two Z^'s cannot vanish at the same time. 

(c) If, for fixed /3, Z/3 = then h2 = h^ = 0. 

Proof. By cyclicity it is sufficient to prove (a) for 
^2 = ^3 = 0, (b) for Z2 = Z3 = and (c) for Zi = 0. 
[&) F2 — F3 = Q forms a linear and homogeneous sys- 
tem in the variables (/12, /13), the determinant of which is 
constantly proportional to x\ and hence cannot vanish. 
Thus /i2 = /13 = 0. 

(b) Elimination of /i from Z2 ^ Z3 = Q yields 2:1 (4a;2 -I- 
7x2X3 + AX3) — 0, contradictory to the Petrov type / as- 
sumption. 

(c) One first calculates that 



02(^1) = (42:^-42:^-|-42;22:3-f/i(2x3-2:2)-MV2)/i2 
+ {4x1 - 42:3 - 42:22:3 + /i(22:2 - 2:3) + /i^/2)/i3. 

Combining the assumptions Zi = and lu°' ^ with 
the equations (|43ll45p it follows that F1F2F3 = 0. 
Now suppose that (/i2,/i3) 7^ (0,0). Then, for fixed k, 
the determinant Dk = Dk (2:2 , 2:3 , /i) of the linear and 
homogeneous system ©2(^1) = Fk — in the variables 
(/i2,/i3) should vanish. As the computed Dk is not a 
multiple of Zi, lemma 13.11 applies with F = Zi and 
G = Dk ■ Because of the hypothesis only the possibility 
jj? — T2 = remains, but elimination of ^ from this 
equation and Zi — yields 2:^2:22:3 = 0, contradictory to 
the Petrov type / assumption. Thus /12 = /is = 0. D 

The key step in the deduction is the following 

Proposition 3.3 UJ1UJ2UJ3 = 0. 

Proof. Suppose on the contrary that LU1UJ21-U3 7^ 0. 
Firstly, if /i2 and /i3 both vanished, one would have 
that aU Fa = and hence, by (Plli5| . aU Z^ = 0, 



which is impossible according to lemma I3.2f b). Thus 
(/i2,/i3) 7^ (0,0), whence also Z1Z2Z3 7^ by lemma 
l3.2r cV Secondly, the equations (|H)) . (|l5l) form a lin- 
ear and homogeneous system in the variables (1^2,^3), 
the determinant of which must vanish by the hypothesis. 
Making analogous observations for the couples (|43|) , (|44)) 
and (I43p . (|^5)l one arrives at 



(48) 



4Za + lZa-ll^a " Fa+lFa^lX^ — 0. 



The key steps are now the following. Substituting 
into the first component of C'-'^'q = one obtains an 
equation of the form 

{Piujj + P2UJI + P3L0I + Qxlxlxl)uji - RX1UJ2LO3 = 0[,49) 

where Pi, P2, P3 are homogeneous polynomials in 
(2:2,2:3,/!) of total degree 7, while Q and R are homo- 
geneous in {x2,X3, fj,) and quadratic, resp. linear homo- 
geneous in (/i2,/i3). Now multiply (|^ with 2Zi, add 
PLi, divide the result by ui and finally eliminate the 



oj^'s by means of p5|) . Performing the analogous oper- 
ations on the second and third components of C^^^^q = 
one derives three equations Ga = 0, where 

Ga = Qa+la+lh2+Qa + la-lh2hz+Qa-la-lh^ + 2ZiZ2Z^^ 

Q22, Q23 and Q33 being homogeneous polynomials in 
(x2,a;3,/i) of total degree 5. Consistency of these equa- 
tions with (/12, /13) 7^ (0, 0) and ZiZ2Z^ 7^ requires that 
F = 0, where 

F = -1107/ -I- T2ii^ + SSTaM^ - 1510r|/ - eSTzTg/i^ 
-h2(130T| - 3T|)/^2 _ UTsTin - 2T2(4T| - T^), 

as is readily deduced by putting ft,2 equal to 1 in G2 — Gi 
and G3 — Gi, and then computing the resultant w.r.t. h^ 
of the resulting polynomials, which yields Z1Z2Z3F — 0. 
By repeating the same procedure for G2 — Gi and Og (F) , 
which is linear homogeneous in (/i2, h^) by p3|) and (|46t - 
I47p . one arrives at a relation G(a;2, 0:3, /i) = 0, where G 
is a homogeneous polynomial of total degree 18 which is 
not a multiple of F. With this F and G one deduces from 
lemma [01 that ^^ — T2 = 0. Calculating the resultant of 
^2 - T2 and i^, however, one finds T^{T^ - T^f = such 
that the Petrov type is O or D, which yields the desired 
contradiction. D 

Proposition 3.4 The vorticity is parallel to a Weyl 
principal vector, i.e., it is an eigenvector of -E^b. 



Proof. By proposition [331 we must show that the case 
wi = 0, 0^21^3 7^ is inconsistent. From this hypothesis 
and (gS]) we get Zi == 0. Hence, by (gl]) with a = 1, 
X2 7^ 2:3 and 



11 = 2- 



X^i - X2 



(50) 



The a = 1-component of (|42p immediately yields (T23 = 0, 
while substitution of ([50)) into the 2- and 3-components 
gives 



(712 = W3- 



X2 + 3X3 
2^2 -Xz 



fl3 = ^2- 



3X2 + X3 
X2-XS ' 



(51) 



respectively. Likewise, the 23-component of ([T4[) yields 
fii = 0, while substitution of ([5D[) into the 31- and 12- 
components gives 



0,2 = 4a)2a;2 



3^2 +x^ 
[X2 - X3)2 



r^s = 4LJ3X3 



2^2 +3:3 
(X2 - X'i)'- 



, (52) 



such that $12^13 7^ 0. Next, by (I44M5P and lemma 
I3.2r a). or by Zi = and lemma [32Ic), we conclude that 
h2 = hz = 0. Propagating this along the dif' integral 
curves, on using the diagonal components of p2[). and 
then substituting (|5ip and ([52| one gets 



2 _ (X2-H2X3)(X2 -X3)^ 

■^2 ^ 48x2(x2 + a;3)^ 



(2x2 +a;3)(x2 - X3)^ 



48X3(X2-|-X3)2 



such that (x2 -I- 2x3) (2x2 + x^) 7^ 0. Applying Oi on the 
left and right hand sides of these two equations one finds 
9uo2 = OloI = 0, whence = 0. With the so far obtained 
equations Raychaudhuri's equation (|lip reduces to 



2x2X3 -I- A(x2 - X3) = 0. 



(54) 



Using ([nni), (IS2D, (ESI) and ([SID one derives the simple 
relations 



^ = 2A + 2(x3-X2), 


(55) 


X2-X3=3(0^+172^A), 


(56) 


X2+X3=i{Ol~Ol). 


(57) 



At this stage ([T2|) -([T5 l) reduces to the vanishing of 
cis, o'i2, 072,^3, X2,X3 and jjL under (?0j such that the 
above algebraic relations either propagate consistently 
along the 9o" integral curves, in a trivial way, or lead 
to (9o^2 — do^z = 0. However, propagating them along 
the da"' integral curves will lead to a contradiction. Do- 
ing this for (|55p and mixing up with equation (J19p and 
the off-diagonal components of ([20| leads to 

9iX2 = 2riX3 + giX2, 9iX3 == -riX3 - 2giX2, (58) 



^2X2 = -^2X1 + q2a;3, ^2X3 = -r2Xi 



'72a;3 



(59) 



^^32^2 = 932:1 



r:iX2 



9^x3 = qsxi - ^3X2. (60) 



The diagonal components of ([^D)) on the other hand re- 
duce to 

(x2 +X3)ni +X2n2 = 0, (x2 + X3)ni + X3n3 = 0. (61) 

Two more algebraic relations follow by suitably combin- 
ing the 2- and 3-components of ([T6)) with the 31- and 12- 
components of ([TT]) . respectively, and substituting ([ST]). 
namely 

W3(x2 + x^)qi + 0)2(3x2 -I- X3)n2 -|- 2^2X271.3 = 0, (62) 
^2(2^2 + 2;3)ri -1- 2013X371,2 + ujz{x2 + 3x3)713 = 0. (63) 

Propagation of ([M]) along the c^i" integral curves leads 
to 



X2(2x2 + x\)qi + X3(x2 + 2x3)ri = 0. 



(64) 

The equations (|6i p -(|64 p form a linear and homogeneous 
system in qi , ri and the 77,^ . Calculating the determinant 
of the system matrix and substituting (|53p one finds (x2 — 
2;3)^(2;2 + 2x3)(2x2 -I- X3)(x2 -l- X3)(x| -|- x§)/16, which 
cannot vanish. Hence 

gi = ri = Til = 7i2 = 773 = 0. (65) 

Next, applying 82 and % to ([54]) leads to 



(53) 



X3(x2 -I- 2x3)172 + 2(x3 - X2)(x2 + X'ifr2 = 0, (66) 
2(x3 - X2)(x2 + xsfqs - X2(2x2 -I- X3)r3 = 0, (67) 



respectively. The comoi(a;2 — x^) commutator relation 
becomes 



^32^392 - t^2X2r3 = 0. 



(68) 



The comoi(a;2+a^3) commutator relation is a combination 
of (|66|) . (|67p and ((68| . However, by propagating the first 
equation of ([55]) along 82'^ and the second one along 83'^ 
we find 

(a;2 - X3)2 
192(a;2 + X3yx2UJ2 
(x3(5a;2 + 2:3) (a;2 + 2x3)^^2 - -^ixj + 4^X2X3 + xl)r2) , 

(x2 - 2:3)2 
192(a;2 + 2:3)3a;3a;3 

(4(2:2 + 42:2.2:3 + 2:3)173 + 2:3(2:2 + 52:3)(2X2 + 2:3)2r3) , 

respectively. Substituting this in (jlSp one gets a new 
independent relation which, together with (|66ll68p . forms 
a linear and homogeneous system in q2, r2, 93 and r^. 
After substitution of (|55|) the determinant of the system 
matrix is a homogeneous polynomial in (x2, 2:3). If it does 
not vanish then (72 = ?'2 == 93 — ''3 = 0; if it does vanish 
then it follows in combination with (|54p that 2:2 and 2:3 
are constant. Looking at (|59| and (|60| we conclude that 
in any case 



92 = »'2 = 93 = T3 = 0. 



(69) 



Finally, substitution of (|5T]|551) . (jSSl) . (P]) . (Pjl and g = 
^2 = ^3 into the 11-component of psp yields 2:2—2:3 = 3 A 
and thus, by comparison with ((56|) . fij + ^3 = 0; which 
is in contradiction with $12^3 7^ 0. D 



By proposition 13.41 the vorticity vector is parallel to 
e.g. di"", such that uJi ^ and 

W2 = CJ3 = 0. (70) 

Then the 2- and 3-components of ([^^ simply read 

a'12 = 0-13 = (71) 

and the 31- and 12-components of (fHl) yield 

O2 = r!3 = 0. (72) 



With these specifications, the equations (|44II45[) are au- 
tomatically satisfied, whereas (H5|) gives Fi — 0, i.e., 



^22:3 + ^132:2 = 0. 



(73) 



Now the 1-components of C'^^-'a and C'"^'o/(12wi) are of 
the form 

fiuji+Fi, fi{h2h3 + als) + F2, (74) 

where Fi and F2 do not contain /i. Calculating one fur- 
ther derivative 03{C^)a and dividing the 1-component 
by 2cji one gets a condition of the form 

{h3~h2Mh2h3 + al-i) + 9ujf~l2{h2X3 + h3X2)]tl + F3 = 0, 



with F4 again independent of fi, which thus can be obvi- 
ously eliminated by means of (|75|) and ([Ti]) . But doing so 
one miraculously arrives at {112X3 -I- h3X2)Fi + X2X3{h3 — 
h2 + 0) = 0, such that 



= ft,9 — ft.3 



I.e. 



11=0. 



(75) 



by the Petrov type / assumption. Thus the restriction 
to u-^ of the expansion tensor O'^i, has a zero eigenvalue, 
9i = ai + 6/3 = 0, as is seen from the representation 
matrix 















~h3 


1723 





^23 


h2 \ 



(76) 



[dafs] 



w.r.t. the WP triad [di" , 82" , 83") . Applying do to ([75 
one gets 



/i = 2(2:2 - 2:3 -I- A) 



A* 



-6E1 + 2A. 



(77) 



The 82- and 93-derivatives hereof, together with 2:22:3 7^ 
0, give 



92 = n3 = 0. 



(78) 



Together with (|75p this implies that the space-times pos- 
sess a geodesic space-like WP vector field, in casu c^i". 
In combination with ([77]) the 23-components of ([H]) and 
(|42)) yield 



CT23 = 



2:2 — 2:3 + 2A 

Xl 



-^1, 



„(22 + A)(x3-A) 

iil — — Z 5 tJl. 



2:? 



(79) 

A further 9o-derivative of ([77]) leads us back to (|73p . 
Propagating ([75)1 along Sq'^ and substituting ([7^ gives 

8AwJ(2;2 -I- A) (2:3 - A) - 22:22:3 (x? - A^^) = 0. (80) 



On using ([73|) and the first equation of ([79|) this may still 
be simplified to 



(/l2^3 + 0'23 - ^l)A + 2:22:3 = 0. 



(81) 



As the Petrov type is assumed to be /, ([50]) or ([5T|) implies 
that A 7^ 0. The 9o-derivative of (|8ip is found to be 
identically satisfied under the already derived equations. 
Equivalently, one may start with the covariant expres- 
sion of (|75p. namely 



7afcW 



i.e. (Tafcw 



-OUn 



with w" 7^ 0. Two covariant time derivatives lead to 

EabLo^ = --(M~2A)cja, 
6 

tr(ai?) = ^(^-2A), 



(82) 
) 

(83) 
(84) 



which, given ^, is nothing but ^^ and dTT]). Now 
C'^^a = simplifies to 



[a,E]a ^-(pL + 2K)L0a. 



(85) 



the covariant time derivative of which may be reduced 
by the aheady obtained equations to 

tiE^ - -in - 2Kf - 2lY{a^E) + -etr{aE) 

D o 

+ — {^i + 4A)(3trcr2 - 26|2) - 2Aw"wa = 0, (86) 
18 

The time propagation of ([55)) onto the WP tetrad is iden- 
ticahy satisfied under the projections of (|82ll86p . 

We conclude that, for perfect fluids satisfying ([21111), 
the complete set of covariant equations generated by time 
propagation of 1121]) is equivalent to i82Vl86]) and, after 
projection onto the Weyl principal tetrad, to ( [70[|7i5[ ), 
^, ^, {^ and ISIP. 



4. FURTHER INVARIANT PROPERTIES AND 
SUBDIVISION 

From here it is advantageous to use 

Ui = r203 — UJl+ (723 , Vl = — r302 ~ LOi— 0-23 (87) 

as auxiliary variables. This allows to rewrite (|79p as 

(a;3-A)t/i = {x2+h)Vi = xiVli = -[(a;2+A)+(a;3-A)]rJi, 

. . (^^) 

which is a linear and homogeneous system in the variables 

{x2 + A, ^3 — A). As we would get xi = —{x2 + X'i) = 

if both these variables were zero, the determinant of the 

system should vanish: 



UiVi + Uini + Vini = Q, i.e. ni = 



-^23 



2^1 



-• (89) 



After some manipulation, the projections of the Bianchi 
identity and Ricci identity for u° reduce to 

doX2 = ~h2{x2 + A), 90^3 = /i3(a;3 - A), (90) 

d^h2 = -hl-2^^Ui-X2, (91) 

5o/i3 = hl + 2J7iVi - x-i, (92) 

00^1 = -Ouji, ^0^23 = -^ihi - ea23, (93) 

9iX2 = -qiX2, dixz = rixz, (94) 

92X2 = -^2X1, dzxz = q3,xi, (95) 

dih2 — -qih2 — 2nicr23 + n^Vi, (96) 

dihs = rih3 — 2nia23 — n2Ui, (97) 
diUi = -qiUi + 2nih2, 91^1 =riVi- 2nih3{98) 

d2h2 + dsVi = -r2hi + 2930-23, (99) 

93/13 + d2Ui == 93/11 + 2r2Cr23, (100) 

qiVi-riUi + 2ni{h2~h3)^0, (101) 

nixi = n2X2 = n3a;3 = 0, (102) 

while the ll-component and subtraction of the 22- and 
33-components of ([35l) yield 



diqi = -q? - 2nin3 - {x2 + A), (103) 

diri=rl + 2nin2-{x3-k). (104) 



On calculating the covariant time derivative of the re- 
maining primary integrability condition (|20p . one checks 
that its projection w.r.t. the WP tetrad is identically 
satisfied under the above equations. For this purpose 
one uses the 'curl-dot' 1-1-3 covariant commutator rela- 
tion applied to the two-tensor Eab^ which in our situation 
becomes j45| 

(curlS')„fc = CWA S ab - IzOcVltX Sab- (J e''^cd(a^''Siy)^ 

-D(a'S'6)''wd -1- W(adiv 5b) . 

Further properties of A may be deduced. Firstly, 
[cr,E]a = (<=> (723 = 0) is impossible: 0-23 — Q ^ tJi, 
([931 and ([89l) would give —hi = /i2 + /13 = 0, and then 
(|91ll92p would yield xi — Q and hence Petrov type D. 
Secondly, the shear tensor, or equivalently the expansion 
tensor, cannot be degenerate. As cr23 ^ 0, this would 
imply ^13/13 -I- (713 = 0, which expresses the vanishing of 
the determinant of the non-trivial 2 by 2 block in ()76p , 
But taking the 9o-derivative hereof and using ([73| , ([89l) 
and X2a;3 ^ would then yield /i2 = /13 = 1723 = 0, a con- 
tradiction. Thirdly, ([77p implies that for each member of 
A the number ig of functionally independent zero-order 
Riemann invariants is at most 2. The following proposi- 
tion states this much more precisely. 

Proposition 4.1 Within A one has 1 < io ^ 2, with 
moreover Iq = 1 '\i and only if one of the three mutually 
exclusive possibilities 6' = 0, t/i = or Vi = is satisfied. 
The cosmological constant is positive for UiVi = u}\ — 
crlg — and negative for — 0. 

Proof. We have t^ <\\i and only if the differentials 
da;2 and dxg are algebraically dependent at each point. 
Then, in particular. 



9oa;2 Sixa - (9ia;2 3^X3 = 0. 



(105) 



From ^, ^, ([MD and a;2a;3 7^ it follows that this 
equation is identically satisfied in the case where 9 = 
/i2 ~ h3 — Q (which implies /12 = /is = 0), and that it 
gives 



(a::3-A)gi+ri(j;2+A) = 



(106) 



when 9 ^Q. In combination with (HH) and xi ^ 0, ([T06l) 
is equivalent to 



giVi+riC/i=0. 



(107) 



Now equations ([TS)) . ([HJ) and the first equation of 
may be solved for X2, X3 and A, giving 



X2 = 



EG 

2/13W1 



X3 



EG . EG"^ 

A 



2h2UJi ' 



4/l2/l3wf 



(108) 



where 



E = /12/13 - UiVi 7^ 0, G = h2Vi + h3Ui ^ 0. 
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From ^ and (I10mi02|) one obtains 

qiVi-riUi qiVi-riUi 



"-2 



2/l2 



"^3 



ni 



2/^3 

qiVi -riUi 
2{h2 - hs) ■ 



(109) 
(110) 



Taking the ai-derivative of (fTOT]) . using JMl) and (fT03l 
HMl) . and substituting pUS)) and (|109mi0p . one surpris- 
ingly finds 



2fi,2n.3t'cjf 



(111) 



Together with ^08]) and 6'a;2X3 7^ it foUows that UiVi = 
0. 

Conversely, if = 0, i.e., /12 = /i3 = 0, then ([9T]l92l) 
implies 2:2 = ~2fliUi and 2:3 = 2riiVi, such that in 
particular UiViili 7^ 0. Inserting this into (jSTjl one finds 
that A = — 4J7^ < 0; whence, without loss of generality, 



A = V^, Qi = 



A 



Ui 



Vi 



X3 



X2 

-' ''^ A 
Herewith, either of the equations (|88p becomes 

2X2X3 + (X2 - X3)A^ = 0, 



(112) 
(113) 



which establishes that the zero-order Riemann invariants 
are algebraically dependent. On the other hand, if e.g. 
14 = ((723 — ^1) then Ui ^ 0, and we consecutively 
deduce X3 = A from ([55)) . X2 — —^2/13 from ([HTjl and 



A = /i§ > from ((73)l . The case f/i = is equivalent to 
Fi = (switching of 2- and 3-axes) and leads to X2 = 
—A, X3 = /i2^3 and A = /i| > 0. 

Finally, if to was allowed to be zero, then h2 — h^ — 
by ([90| . (|73|) and X1X2X3 7^ 0. Taking the combination 
2^3 ^igi - X2diri, using (|103m04p and then (fTU^ . and 
finally inserting gi = ri = (as follows from (|M|) ) we 
arrive at (n^ — A)xi = 0, contradictory to the fact that 
A < when h2 = hj, = Q. D 



It follows from the proof that when to = I7 the 
only allowed functional relations J-'(x2,X3) — are 
X2 = iA'^, X3 = ±A^ or (|113p . with A a real constant. 
Thus, e.g., linear relations involving both X2 and X3 are 
inconsistent. In particular, all members of A are of 
Petrov type /(M+) in the extended Arianrhod-Mclntosh 
Petrov classification [2^|, i.e., the eigenvalues Ea are all 
non-zero. 

For the non-expanding subclass of A one derives from 
([77|. (PI) and (fTT^ that 

w = ~8r2i(r2i + t^i) == 2(tj? + crfg) (i-^j . (114) 

This is positive if and only if |a;i| > |f2i| and ^iloi < 0, 
or if and only if 



where iVab = U[a;b] = ^abc'^'^ IS the vorticity tensor. Thus 
the energy density is positive if and only if the norm of 
the vorticity tensor is larger than the norm of the shear 
tensor. 

On the other hand, the expanding tg ^ 1 family is 
characterized by \uji\ — \a23\- When e.g. Vi = (0-23 = 
LUi) one may take h^ = —X without loss of generality, 
which gives X2 = A/12 and X3 = A^ > 0. Still note from 
()90[) that X3 = const is an equivalent characterization 
for this situation. We further deduce fii = from (|88p . 
n = 93 = from (|94ll95p . m = n2 = 713 = from ^ 
and (|102p . and finally ^2 = by taking the 92-derivative 
of X2 = A/i2 and using (|95p . ([M)) . Herewith the remaining 
equations are identically satisfied or determine consistent 
expressions for derivatives of the remaining variables /i2, 
LUi and Qi. 

The invariant integration for this family was performed 
in [29]. The result is the line element 

X^ds^ = -{dt-2ydzf + {dx + fi{z)dzf 

+e^'[dy+{f2{z)+y^ + e-^')dzf 

+ icosx-e-*ff3{zfdz^, (116) 

where A = vA plays the role of a constant scaling factor 
and where t, x and y are invariantly defined coordinates. 
The arbitrary scalar functions /i(z), f2iz) and 73(2:) are 
invariantly defined; only when these are all constant there 
is a continuous isometry group, which is one-dimensional 
and generated by d/dz. Writing R = cosxe* one has 

A 2A^ / 1 \ 

- = S(I)«' "'TT' "'H' + SJ' <"" 

such that the matter density and expansion scalar are 
positive for i? > 0. 



5. LIMIT CASES 

The special gravito-electric irrotational dust solutions 
found in [l3|, namely 



Ms^ = -dt^ + dx^ + [e 



gi{y)cosxf dy'^ 



[e + g2{z)ii\Tix\ dzJ 



-* '-e2*dz2 



Ms^ = -dr + dx^ + [e 

+.9i(2/)cos(x + 52(y))]^ V 



(118) 
(119) 



|wi|>|ct23| ^ \nJ^^> v^Ibo^, (115) 



were precisely characterized by having a zero eigenvalue 
for the expansion tensor or, equivalently, by possessing 
a geodesic space-like Weyl principal vector field. In the 
above studied rotating case these properties have been 
shown to be satisfied automatically, the singled out WP 
vector field being moreover parallel to the vorticity. 
Thus, surprisingly, A consists exactly of the rotating 
generalizations of illSVlll^) . which are on their turn 
contained as the irrotational limit solutions. The metric 



(I116P is the rotating generalization of (|119p . As also 
follows from the analysis in [ij], the relations (|70ll73p . 
(|75ll78p and (J79ll8ip with wi = are valid for the metrics 
(|118mi9p when taking ^ - 9o and ^ ~ ^i. Note from 
([80| that A = A^ > 0, where A again plays the role of an 
overall scaling factor. 



The shear-free limits (T23 — h2 ^ h^ — oi A are non- 
expanding and automatically satisfy (jllSp . We take ([SS]) 
- instead of ([79]) - as a defining relation for A, and derive 
from (??), dMl), (fTOT|) . (fT02l) . dMl) and ([771), respectively, 



that a;i = 0, n2 = ns = gi = 7'i — 0, X2 — n + A 
-A > and 2fli = -wi 
to the Godel solution 



A. This corresponds precisely 



2Ads^ = dx^ + dj/ 



-e^^dz^ 
2 



(dt-fe^dz2), (120) 



which may indeed be interpreted as dust p — with 
cosmological constant A. 

Nowhere in the reasoning leading to the results of the 
previous section we have explicitly used that /i 7^ 0. Be- 
cause of (|77p and the inconsistency of linear relations in- 
volving both X2 and 0:3 (cf. supra) it follows that Petrov 
type / A-terms cannot be contained in A as 'limits' ji ~ 0. 
Thus we have established: 

Theorem 5.1 Purely electric, algebraically general Ein- 
stein spaces for which the time-like Weyl principal vector 
field is geodesic and rotating do not exist. 

Mars [31| characterized the Kasner space-times [S^, |3j| 
as the purely electric Petrov type / vacua (^ = A = 0) 
with non-rotating and geodesic time- like WP vector field. 
Combining the above theorem with this result and with 
further work on the non-rotating Einstein space case fi — 
7^ A 46], the following conjecture may be stated: 

Conjecture 5.2 The only purely electric, algebraically 
general Einstein spaces with a congruence of freely falling 
Weyl principal observers are the Kasner models and their 
generalizations including a cosmological constant. 



6. CONCLUSION AND DISCUSSION 

The following theorem summarizes the results obtained 
in this paper. 

Theorem 6.1 Consider the class A of algebraically gen- 
eral, gravito-electric rotating dust space-times. For any 
member of A, the vorticity vector at each point is paral- 
lel to a Weyl principal vector which is moreover geodesic, 
the shear tensor does not commute with the Weyl electric 
tensor and is non-degenerate, the Petrov type is /(M+) 
in the extended Arianrhod-Mclntosh Petrov classifica- 
tion, and the cosmological constant cannot vanish. The 
curvature varibles Eab, M, A and kinematic quantities cTat, 
Lu"" and 9 are subject to the covariant algebraic relations 



(I82ll86p . The number to of zero-order Riemann invariants 
is cither 1 or 2. The subclass corresponding to to = 1 
splits into two separate parts. The first part consists of 
all non-expanding solutions [6 — 0); these are Petrov 
type /, shearing generalizations of the Godel universe, 
with A < 0, and the energy density is positive if and 
only if the norm of the vorticity tensor is strictly larger 
than the norm of the shear tensor. The second part cor- 
responds precisely to the solutions for which these norms 
are equal, where now A > 0. The class does not allow for 
Einstein space (A-term) limits /x = 0. 

Note that these statements have been derived for dust 
(j) = 0) space-times with cosmological constant, but they 
remain valid, mutatis mutandis^ for perfect fluids with 
constant pressure. E.g., it has been proved that if an 
algebraically general, gravito-electric, non-expanding but 
rotating perfect fluid model has constant pressure p, then 
p must be larger than the cosmological constant present. 

A second remark concerns the eventual relationship 
between the rotation of a congruence of observers u" in 
a general space-time and the magnetic part of the Weyl 
tensor with respect to it (SJ]. Speculations about such a 
connection stem from the fact that the Bianchi 'div H' 
equation (for perfect fluids or A-tcrms: equation pT|) 
with dbf Ha instead of in the right hand side), contains 
the 'angular momentum density' source term {ji -\- p\jja 
- in contrast to the analogous Maxwell equation [3a, l3g| . 
The link was verified and affirmed, in some sense, for 
e.g. the van Stockum solution [33] [43] and the Bondi 
space-time [38]. However, examples are known of rotat- 
ing gravito-electric perfect fiuids 39, 40]. For such space- 
times the third term in pT|) is exactly balanced by the 
first two terms. E.g., in the case of the Stephani-Barnes 
gravito-electric dust space-times of Petrov type D men- 
tioned in the introduction, the rotation lies in the E°'\j- 
eigenplane and one has 3Eabi^^ = — [c, -E]a = — MWa/2. 
For the metrics of class A considered in this paper, the 
rotation is an eigenvector of E°'i, and there is an exact 
balancing based on ([5^ . (|551) and (|55)) . In this respect, 
/i = — 2A again leads to the homogeneous Godel solution, 
where now the rotation vector of the dust is an i?"h- 
eigenvector corresponding the non-degenerate eigenvalue 
2A/3 and spans the axis of local rotational symmetry. 

Another consequence of the analysis is the follow- 
ing restatement of the generalized silent universe con- 
jecture (as put forward in [14]): an algebraically gen- 
eral gravito-electric dust space-time is either an orthog- 
onally spatially homogeneous Bianchi type I (OSH BI) 
Saunders model flM], or possesses a geodesic space-like 
E'^fj-eigenvector field (and then belongs to A or the non- 
rotating limit families U1^119\) ). Conjecture 15.21 states 
that only OSH BI models are possible for the correspond- 
ing A-term limit. 

By this investigation, the question by van Elst et al., 
whether A constitutes a broad class of solutions (cf. the 
introduction), may be answered affirmatively. We have 
seen that A may be partitioned into three subclasses, 
one constituted by the non-expanding members, one by 
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the expanding members with to — 2 and one by the ex- 
panding members with tg = I. Equation (|116p gives the 
general line element corresponding to the third subclass, 
which provides a first explicit example of algebraically 
general rotating dust with an at most one-dimensional 
isometry group, depending on three free functions of one 
coordinate. 

On the other hand, the non-linearity of some of the 
class- defining algebraic relations hinders a transparent 
consistency analysis for the other two subclasses, at least 
in an orthonormal approach based on an iJ^b-eigentetrad. 
However, as somewhat hidden in such an approach, it 
turns out that there is a significant geometric duality 
between the dust four-velocity m° and the normalized ro- 



tation vector v'^ = w^/yi^-W^ at each space-time point: 
both are geodesic Weyl principal vectors, the vorticity 
vector of the one is parallel to the other [43] , and the one 
is an eigenvector of the shear tensor of the other. These 
properties are most naturally expressed within a 1-1-1-1-2 
covariant formalism, the first '1' standing for m° and the 
second '1' for w", whereas '2' expresses that one leaves 
a 50(2, R) rotational freedom in the orthogonal comple- 
ment of these vectors at each point. There is good hope 
that, on using a complexified version of such a formalism, 
one can elegantly tackle the general consistency problem 
for A, and thereby substantiate that also the remaining 
subclasses to = 2 and to = 1, 9 = contain a large num- 
ber of metrics. 
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A recent attempt to prove this conjecture [1^ contains 
conceptual mistakes, as will be commented on elsewhere. 
To the best of the author's knowledge. 
The Qa are the non-zero components of the angular ve- 
locity vector of the triad {da) w.r.t. the 'inertial compass', 
see e.g. [24j| and references therein. 
See e.g. ^^, pp. 158. 

This commutator relation was given in [23], formula 
(A18), for the further subcase uj°' = 0. In the general 
formula (B.14) of [2y], for imperfect fluids and with Hat 
and u" not necessarily zero, the terms corresponding to 
-\-3Hci^aSb)'^ have the wrong sign. 

In collaboration with N. Van den Bergh it has meanwhile 
been shown that there should be at least an Abelian G2 
isometry group. This will be presented elsewhere. 
The rotating dust interior of the van Stockum cylinder, 
when staying close to the axis, is purely electric w.r.t. a 
certain congruence of non-com,oving observers. 
The vorticity of di°' may vanish. 



